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Abstract 

The explicit formulae for the height distribution and expected number of local max¬ 
ima have been obtained for isotropic Gaussian random fields on certain low-dimensional 
Euclidean space or low-dimensional spheres. 
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1 Introduction 

It has been known that, by using the Kac-Rice formula, one can obtain implicit formulae for 
the height distribution and expected number of local maxima of smooth Gaussian random 
fields parameterized over an A-dimensional space. See, for example, Adler and Taylor (2007) 
and Cheng and Schwartzman (2015). In particular, if the Gaussian field is isotropic, then 
the application of certain Gaussian random matrices involving Gaussian orthogonal ensem¬ 
ble (GOE) makes the above implicit formulae computable [cf. Cheng and Schwartzman (2015) 
and Fyodorov(2004)]. However, for an arbitrary dimension N, the general explicit formulae are 
still hard to formulate due to the difficulty of evaluating certain GOE computation, see The¬ 
orems o and lun below. Motivated by statistical applications of detection of peaks in Cheng 
and Schwartzman (2014), we investigate here the explicit formulae for the height distribution 
and expected number of local maxima of isotropic Gaussian random fields on A-dimensional 
Euclidean space or A-dimensional spheres, where A = 1, 2, 3. 

2 Preliminary GOE Computation 

Recall that an A x A random matrix Mjv is said to have the Gaussian Orthogonal Ensemble 
(GOE) distribution if it is symmetric, with centered Gaussian entries Mij satisfying = 

‘Research partially supported by NIH grant R01-CA157528. 


1 



1, Var(Mjj) = 1/2 if i < j and the random variables {Mij, 1 < i < j <A^} are independent. 
Moreover, the explicit formnla for the distribntion Qn oi the eigenvalues Xi of Mjy is given by 


N 


Q7v(dA) —— JJe ^^idXi — Aj|1{Ai<...<Ajv}5 

l<i<j<N 


CN , 
2 = 1 


( 2 . 1 ) 


where the normalization constant cjv can be computed from Selberg’s integral 

N 

■ 2 , 


c« = i( 2 ^ 2 )«nr(i + | 
2=1 


( 2 . 2 ) 


Denote by AA(/i, cj^) a Gaussian distribution with mean /x and variance cr^. Let <l>(x) = P(X < 
x) and <hs(xi,X 2 ) = P(Ali < xi,X 2 < X 2 ), where X ~ AA(0,1) and (Xi,X 2 ) is a centered 
bivariate Gaussian vector with covariance S. For 7 > 0, let 


Gy 





(2.3) 


It can be seen that Gy{x) = Y^ 7 r/ 7 <h(v/ 27 x). 


Lemma 2.1 For constants 7 > 0 , (T>0, Q!>0 and /3 € M, 


n e-^^^-^^^Gy{x)dx = , 

J-00 Vv« + 7 / 

/ CO py TT^I'^ 

^-^{y-0? / e-'^"'G.,(x)dxd2/ = ^=<I>s(0,/3), 

-CX3 J —00 


where 


{Zi - Z 3 ) ~ Af{-l3, ^) and hence 


27 ^ 


/ g-+7 

_J_ \ 

— 1 20-7 

2 ct I 

1_L 

tr+a / 

\ 2t7 

2fTQ / 

27) 

and Z 3 


(2.4) 


2 a > 


/ CO 

g{x)e-^^^-^^"dx = ^P(^i - Z 3 < 0 ) = ^=<^> 

-00 


vr / l3y/2^\ 


y/oPy ' " " ' V\/« + 7 y ' 

Let g = (0, —/I) and let S be as defined in (|2.4I1 . Then (Zi — Z 2 , Z 2 — Z^) ~ AA(/x, E), which 
implies 

vr3/2 


/ CO py 

g-a(y-/3)2 / g(^x)dxdy = -^-=¥iZi - Z2 < 0, Z2 - Z3 < 0) 

-00 J-oo 

"'^':‘l>s(0,/3). 




□ 
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Lemma 2.2 Let N = 1. Then for constants a > 0 and 6 E R, 


E 


N+l 

GOE 



:^N+i - o.{Xn+i - 



\/ 4a + 2 ab^ 

- -e 2“+i + 

4a 



\/2a + 1 y 


Proof By (12.111 and integration by parts, 


CAr+iE^+^ <j exp 


2^‘n+i ~ O'i^N+i — by 


} /*00 p \2 

= / e-^(A2-Ai)dAi 

J —OD J —OO 


= b r Gi/2(A2)e-“(^2-fe)^^^2 + + 


2a 


Applying Lemma l2.ll and noting that cat+i = 2^/^T when = 1, we obtain the desired result. 

□ 


Lemma 2.3 Let N = 2. Then for constants a > 0 and 6 E R, 


^GOE i 


o^N+i - o,{Xn+i - by 


= ( - + 26^ - 1 ) _ 

. Q / \/2a \ \/a + 1 


6-v/^ \ 6\/fl + 1 - 


+ 




e “+i 


+ 


g 2a + l ^ 


V 2 ab 


\/2a A 1 


a/ ( 2a + l)(a + 1) 


Proof By (12.ip and integration by parts, 

'1 


CN+1^G0E 1 exp 


-A^+i - a(AAr+i - by 


^-a{X 3 -b)^ 


dXs 


rXs 


a 2 

e 2 


/ 


— OO 
OO 


(A 3 — \2)d\2 f 

J —OO 


e 2 (A2 — Ai)(A 3 — Ai)dAi 




(2A3 — 1 )Gi(A 3) + e 2 (7^/2(-^3) + ^36 


dA3. 


Applying integration by parts and Lemma I^TTl and noting that cn+i = when N = 2, we 
obtain the desired result. □ 
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Lemma 2.4 Let N = 3. Then for constants a > 0 and 6 E R, 


"I 


T-^AT+l “ O'i^N+l — b)" 


24a^ + 12a^ + 6 a + 1 2 6 a^ + 3a + 2 3 

- b + + - 


+ 


2 a( 2 a + l)2 

-b^ + 


4a2(2a + l) 2j y^2(2a + 1) 


ab^ 


-g 2a+l$ 


a + 1,2 . 1 — a 


2a “ ' 2 a 2 

+ ( 6 a + 1 + 


- 1 


1 

' V2ab \ 

\/ 2 (a + 1 ) 

l)(2a + 3) y 


+ 

where 


2 3(1 — a 


28a^ + 12 a + 3\ b 

2 a( 2 a + 1 ) ) 2V^{2a + l)V 2 a + 3 

'/^b 


g 2a+3 




[<I)si( 0 , b) + 4>s2(0, b)] 


2y/2ab 


\J (2a + l)(2a + 3) 



( i 

-i A 


Sr=| 

1 2 

14 

l+a 1 
2a / 

M 

to 

II 


3 

2 

-1 


-1 

l+ 2 a 

2a 


Proof By (12.ip and integration by parts, 


^N+I^qOE i 


o-^Af+1 “ O'i^N+i - by 


/ CO /•A 4 ^2 p 

e-aix,-b) 

-CO J —00 J — 


f —CO 

y\2 


r-As 

e 2 (A 4 — A 2 )(A 3 — X2)dX2 

' —CO 


e 2 (A 2 — Ai)(A 3 — Ai)(A 4 — Xi)dXi 


1 2a + 1 + 4a^6^ b 3b 2ab 3 — 3a 3 
, ^ ~ ( 2 a + 1)2 ^4^ ^^ 2 a + 1 4a2 2 




_^ 

X e 2 a+i / e 

J —CO 


2a-|-l / -x 2at> \2 

2 2n, + l ^ 


+ 

+ 


/ 


2“+!'' G*4(A4)(iA4 

A4 


,-a(A4-fe)2 


e" 2 Gi(A3)dA3dA4 


1,3 3, 

- 'rb'^ - -b 

2a 2 _ 

,1 1 , a + 1 + 2a^6^ , b 3b ab 3 — 3a 

2 2a? 2 (a + 1 )^ 2 a^ 2 a a + 1 4a^ 


X e-f^ r e-(“+i)(A4-f^) 

J 


ab \2 

“+l Gi/2('^4)<^A4 


+ 

+ 


3& 1,3 3^ 

- \-b-^ - b 

2a 2 


/ CO /*A4 

g-a(A4-fe)2 / e-^iGi/2(A3)dA3dA4 

-CO J —CO 


f 


( 6 a +l)(2a + 3) ^ (28a2 + 12a + 3)6 

■A 4 + 


4a(2a + 1) 


4a(2a + 1)2 


3\i 


,-a(A4-6)%-^ 


e 2 dA 4 . 


Applying integration by parts and Lemma l2.ll and noting that cat+i = 27r when N = 3, we 
obtain the desired result. □ 
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3 Isotropic Gaussian Random Fields on Euclidean Space 


In this section, we consider smooth isotropic Gaussian fields on Euclidean space. Let {/(t) : t G 
be a real-valued, centered, unit-variance isotropic Gaussian field. Due to isotropy, we 
can write the covariance function of the field as K{f{t)f{s)} = p{\\t — s|p) for an appropriate 
function p(-) : [0, oo) —)• M, and denote 

P' = P'(0), p" = /(O), K = -p' I (3.1) 

Let fi{t) = and fij{t) = ■ Denote by V/(t) and V^/(t) the column vector 

(/i(t),..., and the N x N matrix (/jj(t))ij=i,...,7v, respectively. By isotropy again, the 

covariance of (/(t), V/(t), V^/(t)) only depends on p' and p”. In particular, Var(/i(t)) = —2p' 
and Var(/jj(f)) = 12p" for any i G {1, • • • ,N}, which implies p' < 0 and p" > 0 and hence 

K > 0. 

Over the unit open cube (0,1)'^ in , define respectively the expected number of local 
maxima of / and the expected number of local maxima of / exceeding level u as 

M(/) = # {t G (0,1)"^ : V/(f) = 0, index(v2/(t)) = N} , 

Muif) = #{te (0,1)"^ : fit) > u, Vfit) = 0, index(v2/(t)) = N} , 

where the index of a matrix is regarded as the number of its negative eigenvalues. Define the 
height distribution of a local maximum of / at some point, say the origin, as 


Fiu) := liniP|/(0) > ri|3 a local maximum of /(t) in (—. 
Then we have the following result in Gheng and Schwartzman (2015). 


Theorem 3.1 Let {fit) : t G M'^} be a centered, unit-variance, isotropic Gaussian random 
field satisfying the conditions (Cl), (C2) and (C3) in Cheng and Schwartzman (2015). Then 
for each u G M, 


/ 2 \ {■^+l )/2 

E{M(/)} = T 


AT-h 1 



N/2 


E: 


N+l 

GOE 




and 


Fiu) 


nMuif)} 

E{M(/)} 

(1 - J” (/)(x)E^+^ {exp [X%+i/2 - (Ajv+i - kx/\/2)V(1 - n^)] } dx 

^GOE } 


where (fix) is the density of standard Gaussian variable and p', p” and k are defined in ()3.1I1 . 
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Remark 3.2 Recall the condition (C3) in Cheng and Schwartzman (2015) that k < 1. We 
show below that, the factors (1 — and (1 — in (|3.2|) will be cancelled out, implying 

that (13.2p also holds when k = 1. In fact, by (j2.in . 

(1 - {exp {\n+i-K x/Vlf/{I-K^) | 

=- n -2 / e JJe 2 \Xi - Xj\ 

CN+lVl - fJi l<i<i<Ar 

X n - XN+l\l{\i<...<\,^<Xr^_^_i}dXi-■ ■ dXNdXN+1 


l<k<N 
1 


CN+1 


N 


[ e ^^+i]Te 2 ^idXi TT \Xi-Xj\ TT |Afc - -x/l - k‘^Xn +1 - kxIV2 \ 
r=i i<ij<iv 


l<A:<Af 


^ ^{Ai<...<Ajv<\/WK2Aiv+i+K3;/V2}'^'^l ' ' ' dXNdXN+1, 


where we have made change of variables Aw+i = (Aw+i — Kxjy/2)ly/l — k?. 

From now on, we denote by h the density of the height distribution of local maxima of /, 
i.e. h{x) = —F'{x) or F{x) = h{t)dt. Also notice that, due to the formula E{M„(/)} = 
F(ri)E{M(/)} in (|3.2p . we only need to find F{u) and E{M(/)}. 


Proposition 3.3 Let the assumptions in Theorem \3.1\ hold. If N = 1, then 


E{M(/)} = 


27r V /o' 


and 


, , , a/3 — —7^^^ 2KXy/TT ,, , / nx 

h{x) = ^ e 2(5^ +- 


a/^ 


V6 




If N = 2, then 


and 


E{M(/)} = - 


y/Sirp' 


h{x) = a/3k^(x^ — l)(/)(a:)<h 


KX 


KX 


V 2 ^J 2 tt 


+ 


V 6 


Y^7r(3 - k2 ) 


g 2(3-«;2)^ 


KX 


a/(3 - k2)(2 - k2) 


If N = 3, then 


E{M(/)} = 


29^6-36/^ p"\3/2 
367r2 
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and 


h{x) = 


144(/)(x) 
29^6 - 36 


K? [(1 — + 6(1 — + 12(1 — K?) + 24] 

4(3 - k2)2 


where 


2(1 - ^2)3 3(1 _ ^2)2 _ ^2) 3 

^ 4(3 - k2 ) + 2 


g 2(3-re2) 


= $ 


2kx 


^k‘^{2 — K^) 2 K^(l — K^) 

—-- - - 1 

4 2 


+ 


v/2(3 - k 2) \^y/(3-K2)(5-3K2)^ 

2 

g~ 2(2 -k^) 


=4- 


KX 


Y ^ 2(2 - «; 2 ) \ y /(2 - k 2)(5 - 3 ^ 2 ) 


+ 


7 _ ^2 ^ (1-_K^(1 - ^2)2 + 12(1 - k2) + 28] 


2(3 - k2) 


KXe 2(5-3 k^) 


4y'^(3 — k2^^5 _ 3^2 




H-^—x(x^ - 3) <i>si(0, kx/\/ 2) + <&S2(0)'^a:/\/2) 


E. . 7^ 


-1 

3— 


, 512 = 


3 _ 1 

2 2 

1 2 -k^ 

2 2 


Proof Let = 1. By Theorem 13.11 we see that E{M(/)} can be obtained by applying 
Lemma 12.21 with o = 1 and 6 = 0, while h{x) can be obtained by applying Lemma 12.21 with 
a = (1 — k2)-i and 6 = KxjV2. The cases for N = 2 and N = 3 can be proved similarly by 
applying Lemmas 12.31 and 12.41 respectively. □ 



Figure 1: Density functions of the height distribution of local maxima on Euclidean space. All 
three cases N = 1 (black), N = 2 (blue) and N = 3 (red) are displayed. Solid lines represent 
K = 1, dashed lines represent k = 0.5 and dotted lines represent k = 0.1. 
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Remark 3.4 The regularity condition (C2) in Cheng and Schwartzman (2015) implies that 
< 3 when N = 1, < 2 when N = 2 and < 5/3 when N = 3. Although Theorem 13.11 

holds for K < 1 as shown in Cheng and Schwartzman (2015), we conjecture that the formulae 
obtained in Proposition 13.31 are valid for any isotropic Gaussian fields satisfying conditions 
(Cl), (C2) and (C3) in Cheng and Schwartzman (2015). 


4 Isotropic Gaussian Random Fields on Sphere 


In this section, we denote by the A^-dimensional unit sphere and let / = {f{t) : t € be a 
real-valued, centered, unit-variance isotropic Gaussian random field on . Due to isotropy, 
we may write the covariance function of / as C{{t, s)), t, s G §^, where C'(-) : [—1,1] —>■ M and 
(•, •) is the inner product in See Cheng and Schwartzman (2015) and Cheng and Xiao 

(2014) for more results on the covariance function of an isotropic Gaussian field on We 
define 

C = C'(l), C" = C"(l), Ki = C'/C”, K 2 = C'VC". (4.1) 

It is known from Cheng and Schwartzman (2015) that C > 0 and C" > 0. 

Let be a geodesic (open) ball on with radius 6 centered at to G . Let do 

be the number such that the area of B{to,So) equals 1. Following the notation in Cheng 
and Schwartzman (2015), over the unit geodesic ball B(to, 6 o) in §^, define respectively the 
expected number of local maxima of / and the expected number of local maxima of / exceeding 
level u as 

M(/) = # {t G B{to, do) : V/(t) = 0, index(V2/(t)) = N} , 

Muif) = #{te B{to, do) : fit) > u, V/(t) = 0, index(VV(t)) = N} . 

Define the height distribution of a local maximum of / at some point, say to, as 


F(u) := limP{/(to) > tild a local maximum of /(t) in i?(to,e)} ■ 
Then we have the following result in Cheng and Schwartzman (2015). 


Theorem 4.1 Let {/(t) : t G he a centered, unit-variance, isotropic Gaussian random 
field satisfying the conditions (Cl'), (C2') and (C3') in Cheng and Schwartzman (2015). Then 
for each u G M, 


E{Mif)} 


V 2 J 


E 


AT-1-1 
GOE 




l-\- Ki 


and 

Fiu) 


nMuif)} 

E{M(/)} 


/ 1 + Kl / 


exp 

-^^+1 i>'N+l—\^x/V^) 

' 

> dx 

1 + Kl—K2 / 

2 

1 + K.I—K2 

^GOE- 

exp 

rx2 \2 1 

^N+1 ^iV+1 

2 1-l-rei 

' 
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(4.2) 




















where ki and K 2 are defined in (HI]). 


Remark 4.2 Recall the condition (C3') in Cheng and Schwartzman (2015) that K 2 — < 1. 

Similarly to Remark 13.21 it can be shown that, the factors (1 + ki — and (1 + aci — ^ 2)”^ 

in (j4.2n will be cancelled out, implying that (14.21) also holds when ^2 — = 1. 



Figure 2: Density functions of the height distribution of local maxima on sphere. All three 
cases N = 1 (black), N = 2 (blue) and N = 3 (red) are displayed. Solid lines represent m = 1 
and K 2 = 2, dashed lines represent ki = 1 and K2 = 1 and dotted lines represent ki = 0.1 and 
K2 = 0.1. 


Now, we can obtain the following explicit formulae for E{M(/)} and the height density h, 
which imply F and 'E{Mu{f)}. 


Proposition 4.3 Let the assumptions in Theorem \4-l\ hold. If N = 1, then 


E{M(/)} = 


2TTy/Kl 


and 


h{x) = 


'\/3 + 


Kl 


--e 2(3+/1 i-k2) -\- y/2/KK2X(p{x)^ 

V 27r 




^/3 + Kl - K2 


If N = 2, then 


E{M(/)} = 


+ 


dvr 27rKi-y/3 + ki 
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and 


h{x) = 


2\/3 + Ki 
2 + ki\/3 + K\ 


[ki + K 2 {x^ - 1 )] 


y/i^X 


^/2 + Kl- K 2 


Y^K2(2 + Kl- K 2 ) - ' 

_|-Y---- -xe 2 ( 2 +ki-K 2 ) 


+ 


2 tt 

V 2 


+ Kl- K 2 ) 


(3+/tj^)a;^ 

-g 2(3+k;i-K2) ^ 


y/K^x 


If N = 3, then 


E{M(/)} = 


7r2^3/2 \2V3T^ 


V^(2 + Ki - K2)(3 + ki - K 2 ) J J 

(ki - 1)(ki + 2 ) 


3 (l + Ki)(2Kf + 7Ki + ll) 
2^ 4(3+ Ki) 


+ 


4.^2 + 


Ki 


and 


h{x) = 




3 (1 + ki)(2k2 +7^1 + 11) 

2^ 4(3+ Ki) 


+ 


(ki - l)(ft;i + 2 ) 

4^2(2 + Aci) 


-1 


\2v^2(3 + ki) 

^2 [(1 + — ^ 2 )^ 4“ 6(1 + Ki — ^ 2 )^ 4“ 12(1 4“ Ki — K 2 ) 4“ 24j 

4(3 4- Ki - ^ 2 )^ 

2(1 +Ki - ^ 2 )^ 4-3(1 4-Ki - ^ 2 )^ 4-6(1 4-Ki - K 2 ) 3 

^ 4(3 + ki -K 2 ) ^2 




V^2(3 + Ki — K 2 ) 




X g 2(3 +ki-K2) $ 


2 y/KfX 


^{3 + Ki - K2)(5 + 3ki - 3K2) 


+ 


K2(2 + Ki-K 2) 2 , - «^2)(1 + Kl - '^2) T 

4 ^ + 2 


- \/ 2(2 + Ki — K 2 ) 


X g 2(2+K]^—^^ 2 ) ^ ^ 

+ 7 Ki — K 2 


y/K^X 


Y^(2 + Ki - K2)(5 + 3 ki - 3 k 2 ) j 

3(1 + Ki — ^ 2 )^ + 12(1 + Ki — ^ 2 )^ + 28(1 + Ki — K 2 ) 




4y'i^(3 + Ki — K2)\/5 + 3 ki — 3^2 


2(3 + Ki - K 2 ) 

Q 2 

OK2X 

2(5+3kj^— 3 ^ 2 ) 


+ [k2 X^ + 3(ki - K 2 )] 


1 °’ V 2 y 




where 


El = ' 2 


-1 


_i 3 +ki—K 2 
2 


, E 2 = 


3 _1 

2 2 
_ ]_ 2 +ki—K 2 

2 2 


Proof Let N = 1. By Theorem 14.11 we see that E{M(/)} can be obtained by applying 
Lemma [2.21 with a = (1 + ki)“^ and 6 = 0, while h{x) can be obtained by applying Lemma 12.21 
with a = (1 + Ki — ^ 2 )”^ and 6 = ^J^xl\f2. The cases for = 2 and N = 3 can be proved 
similarly by applying Lemmas 12.31 and 12.41 respectively. □ 
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Remark 4.4 The regularity condition (C2') in Cheng and Schwartzman (2015) implies that 
K 2 — K,i <3 when N = \, K 2 — <2 when N = 2 and K 2 — ki < 5/3 when N = 3. 

Although Theorem 14.11 holds for K 2 — < 1 as shown in Cheng and Schwartzman (2015), we 

conjecture that the formulae obtained in Proposition 14.31 are valid for any isotropic Gaussian 
fields satisfying conditions (Cl'), (C2') and (C3') in Cheng and Schwartzman (2015). 

References 

[1] Adler, R. J. (1981), The Geometry of Random Fields. Wiley, New York. 

[2] Adler, R. J. and Taylor, J. E. (2007), Random Fields and Geometry. Springer, New York. 

[3] Aza'is, J.-M. and Wschebor, M. (2008), A general expression for the distribution of the maximum 
of a Gaussian field and the approximation of the tail. Stock. Process. AppL, 118, 1190-1218. 

[4] Aza'is, J.-M. and Wschebor, M. (2010), Erratum to: A general expression for the distribution of 
the maximum of a Gaussian field and the approximation of the tail [Stochastic Process. Appl. 
118 (7) (2008) 1190-1218]. Stock. Process. Appl, 120, 2100-2101. 

[5] Cheng, D. and Xiao, Y. (2014), Excursion probability of Gaussian random fields on sphere. 
Bernoulli, to appear, http://arxiv.org/abs/1401.5498 

[6] Cheng, D. and Schwartzman, A. (2015). Distribution of the height of local maxima of Gaussian 
random fields. Extremes, to appear, http://arxiv.org/abs/1307.5863 

[7] Cheng, D. and Schwartzman, A. (2014), Multiple testing of local maxima for detection of peaks 
in random fields. arXiv:1405.1400 

[8] Cramer, H. and Leadbetter, M. R. (1967), Stationary and Related Stochastic Processes: Sample 
Function Properties and Their Applications. Wiley, New York. 

[9] Fyodorov, Y. V. (2004), Complexity of random energy landscapes, glass transition, and absolute 
value of the spectral determinant of random matrices. Phys. Rev. Lett., 92, 240601. 

[10] Schoenberg, 1. J. (1942), Positive definite functions on spheres. Duke Math. J., 9, 96-108. 

[11] Schwartzman, A., Gavrilov, Y. and Adler, R. J. (2011), Multiple testing of local maxima for 
detection of peaks in ID. Ann. Statist., 39, 3290-3319. 


Dan Cheng: Department of Statistics, North Carolina State University, 2311 Stinson 
Drive, Campus Box 8203, Raleigh, NC 27695, U.S.A. 

E-mail: dcheng 2 @ncsu.edu 


Armin Schwartzman: Department of Statistics, North Carolina State University, 2311 
Stinson Drive, Campus Box 8203, Raleigh, NC 27695, U.S.A. 

E-mail: aschwar@ncsu.edu 

URL: http://www4.stat.ncsu.edu/'schwartzman/ 


11 




